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ABSTRACT 

We show that candidate contact binary asteroids can be efficiently identified from sparsely sampled 
photometry taken at phase angles a > 60°. At high phase angle, close/contact binary systems 
produce distinctive lightcurves that spend most of the time at maximum or minimum (typically 
>lmag apart) brightness with relatively fast transitions between the two. This means that a few 
(~5) sparse observations will suffice to measure the large range of variation and identify candidate 
contact binary systems. This finding can be used in the context of all-sky surveys to constrain the 
fraction of contact binary near-Earth objects. High phase angle lightcurve data can also reveal the 
absolute sense of the spin. 

Subject headings: minor planets, asteroids — solar system: general — techniques: photometric 



1. INTRODUCTION 

For most purposes, small solar system bodies are best 
observed at solar opposition. Measurements taken at low 
phase angle (the Sun-target-observer angle, denoted a; 
see Fig. [T]) benefit from strong backscattering, thus ren- 
dering the target brighter. In astrometric or surface- 
averaged measurements, the enhanced signal-to- noise 
due to the low phase angle geometry is desirable. How- 
ever, if the goal is to derive information on the shape 
of the target, then high-a data become useful. This has 
been noted in the past, in the context of lightc urve in- 
version problems (|Durech fc KaasalainenI l2003l ) . As a 
body rotates, the amount of sunlight it refiects to the 
observer depends on the time varying projected cross- 
section, which in turn depends on its shape. If well- 
sampled lightcurve information is obtained at multiple 
observing geometries, then the shape of the object can 
be derived. At low phase angles the non-convex features 
of an object are not apparent from its lightcurve. This 
means that in the absence of high-a data the non-convex 
figure of a contact binary can be wrongly interpreted as 
being composed of a single convex component. 

One known effect of high-a observations on lightcurves 
is the increase in the peak-to-peak range of variation 
(fCehrelsl Il956f ). For single ellipsoidal bodies, and for 
phase angles a < 40°, the lightcurve rang e Am has been 
shown to vary roughly linearly with a (|Zappala et al.l 
Il990f ). This is generally known as the amplitude-phase 
relationship (APR). When observing from Earth, the 
maximum phase angle at which a small solar system body 
can be imaged depends essentially on its heliocentric dis- 
tance. The distant Kuiper belt objects can be observed 
up to a ~ 2°, Jovian Trojan asteroids at a < 10°, and 
main-belt asteroids at a < 30°. Only near-Earth objects 
(NEOs) can be observed at phase angles a > 30° which 
makes them the best candidates for high-a studies. 

On the order of twenty binary N EOs have been discov- 
ered, both photometrically (e.g.. iPravec fc Hahnlll997l : 
|Pravec et al. 199 ^) and from radar observations (e.g., 
[Nolan et all I200QI : lOstro et al.l l2002f ) . From the known 
sampl e, it is possible to note a few characteristic fea- 
tures (jPravec et al.ll2006[ ). The primaries have mean di- 
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ameters of ~ 1 km (with a spread of one order of magni- 
tude), usually two to four times those of their secondary 
companions. One case, 69230 He rmes, has a size ratio 
close to unity (i Margot et al.ll20Q3l ). The separation dis- 
tances are small, usually a few primary radii, and the 
mutual orbits are nearly circular. Most measured pri- 
maries spin (around the shortest physical axis) close to 
the critical break-up rate for strengthless bodies with 
mean density ~ 2 gcm~^. The secondaries have synchro- 
nized spin and orbital periods. In the few (~10) cases 
where shapes can be inferred, the primaries are well de- 
scribed by oblate spheroids, while the secondaries have 
elongated shapes along the line of centers. One NEO, 
4769 Castali a, is suspec ted from radar data to be a con- 
tact binary (Ostro et al. 1990). 

In recent years, explanations to how binary NEOs may 
form have mostly converged on the idea of the splitting of 
a parent body. iWeidenschillingi (jl980) proposed that bi- 
naries may form by rota tional fission o f an ob ject effected 
by an off-center impact. "Margot et alj (|2002l ) argued that 
spin-up due to tidal interactions during close encounters 
with planets was the most likely formation process for 
NEO binaries, and Bottke et al.l (|2006l ) were the first to 
suggest that radiaton forces (YORP effect, Yarkovsky- 
O'Keefe-Radzievskii-Paddack) provide the torque needed 
for rotational fission. The main features of the NEO bi- 
nary sample described above have led Cuk (2007) to pro- 
pose the following formation mechanism: due to YORP, 
~ 1 km objects are spun up close to break-up rates, which 
leads to the accumulation of regolith around the equa- 
tor. As YORP continues to spin up the body, the surface 
material is eventually stored in orbit around the (now) 
primary, where it coagulates into a secondary; the lat- 
ter grows into an elongated, nearly Roche equilibrium 
shape balancing gravitational and inertial accelerations. 
The subsequent evolution of the mutual orbit is driven 
by BYORP, the binary version of YORP, which may lead 
the components away from each other but may also bring 
them together to mutual contact (ICuk fc BurnsI l2005l : 
"Scheeres 2007). 

To obtain a well-sampled lightcurve of an NEO, ca- 
pable of unequivocally establishing its nature, would re- 
quire extensive observations at multiple geometries. For 
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Fig. 1. — Angles used to describe the observing geometry. The 
orbital spin axis is denoted s. The contact binary is assumed lo- 
cated at the vertex where the angles (aspect, the angle between 
the orbital spin axis and the line-of-sight), a (phase, the angular 
distance between the Sun and the observer, as seen from the ob- 
ject), and (f) (azimuth, the sky-plane angle between the projections 
of the orbital spin axis and the object-Sun vector) are measured. 



TABLE 1 

Four close/contact binary systems considered. 







(B/A,C/A)^ 


(6/a, c/a)^ 




H.E.^ 


1 


1.00 


(0.66, 0.60) 


(0.66, 0.60) 


1.00 


yes 


2 


0.67 


(0.77, 0.69) 


(0.53, 0.49) 


1.00 


yes 


3 


0.25 


(0.92, 0.83) 


(0.51, 0.48) 


1.19 


yes 


4 


0.13 


(1.00, 0.90) 


(0.51, 0.48) 


1.00 


no 



^ System number^ Mass ratio*^ Primary and secondary axis 
ratios^ Orbital distance in units of A + Components in hydro- 
static equilibrium 



a sparse survey such as will be undertaken by Pan- 
STARRS (typically revisiting each target ~ 4 times per 
lunation), it may take years before the shape of the object 
can be determined. The method presented here can sin- 
gle out potential contact binaries after ~ 5 measurements 
at high-Qf. In this letter, we present the first light curve 
simulations of contact binaries at high phase angle, and 
use them to address the detectability of such systems. 

2. HIGH PHASE ANGLE LIGHTCURVES 

The characteristic morphological features of contact 
binary lightcurves are the large range of variation, typ- 
ically 0.9 mag or more, the rounded, inverted- U shaped 
maxima, and the sharp V shaped minima (^Zappala et al.l 
1 19801 : iLeone et al.l 11983 ). However, these traits appear 
only at phase angles close to a = 0°. Figure [2] shows 
lightcurves of a symmetric contact Roche binary (system 
#1 from Table 1) at four different phase angles, a = 1, 
30, 60, and 90°. The surface scattering of light is mod- 
elled using a Lommel-Seeliger function, taken to repre- 
sent a low albedo, lunar- type surface. The procedure 
used for generating the model lightcurves of Roche bina- 
ries, as well as of Jacobi triaxial ellipsoids, is described in 
detail in Lacerda fc Jewitt (2007i) . The most noticeable 
mutations as a increases include: 

1. the shape of the minimum changes from V-shaped 




J I I I I I I I I I I I I I I I I I L 



0,0 0,2 0,4 0,6 0,8 1,0 
rotational phase 

Fig. 2. — Change in lightcurve morphology with increasing phase 
angle. Lines are labelled with the corresponding phase angles 
(a = 1°,30°,60°, and 90°). Insets show rendering of the binary 
used to produce the lightcurves (phase angle increasing from top 
to bottom) at rotational phases (j) = 0.1, 0.3, 0.5, 0.7, and 0.9. 

first to flat and then to slightly W-shaped; 

2. the transition between low and high brightness be- 
comes sharper; 

3. the positions (in rotational phase) of the minima 
and maxima drift to the left; 

4. the overall brightness decreases, as the illuminated 
fraction of the surface diminishes; 

5. the trough-to-peak range of variation increases. 
Below, we discuss some of these points in more detail. 

2.1. Lightcurve Shape and Detection Probability 

As a consequence of points 1 and 2 above, the prob- 
ability of detecting the total range of brightness varia- 
tion from only a few measurements increases significantly 
with a. This is shown in Fig. [3] where we plot the prob- 
ability of measuring the extent of variation of a contact 
binary lightcurve for different values of a. The Figure 
was generated as follows: each of the four lightcurves in 
Fig. [2] (corresponding to each a) was sampled five times 
at random rotational phases, and the maximum range 
Ami within the set of five measurements was registered. 
This procedure was then repeated for z = 1 to 10000, and 
the cumulative distribution of ranges Ami was plotted 
for each a, starting from the maximum range. The plot 
shows that beyond a ~ 60° most five-point samples of 
the lightcurve are able to identify a variation larger than 
1 mag. For example, at a = 60° there is ~80% chance of 
detecting a Am > 1 mag from five sparse observations of 
the contact binary considered. 

Figure [3] was obtained assuming the most favorable 
observing geometry (Fig. [1]): 6 = 90°, measured be- 
tween the spin pole and the line-of-sight, and (j) = 90°, 
measured in the plane of the sky between the spin pole 
and the object-Sun vector. The total range of a contact 
binary lightcurve decreases with the aspect angle at 
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about 0.03 to 0.04mag/° (jLacerda fc Jewittll2007t ): the 
maximum range is obtained at ^ = 90°. The azimuthal 
angle (j) also influences the light curve range and shape. 
As (j) moves away from 90° (the direction does not mat- 
ter) the lightcurve becomes asymmetric and the range 
of variation is slightly decreased (Fig. [4^). Therefore, 
requiring a favorable geometry, i.e., angles larger than 
chosen minimum values > 6^ and ^ > reduces 
the probability of detection of large Am. Assuming ran- 
domly oriented spin axes, the probability of having si- 
multaneously favorable and is given by 



47r 



sin((9)d0di9 + 

(tT - ^m) 



(1) 



sin((9) #di9 



COs(^n 



which for minimum aspect and azimuth 0^ = (j)^ = 75° 
is ~ 0.15. This probability must be taken into account 
when estimating the intrinsic fraction of contact binary 
systems from a sparsely sampled survey. At any rate, 
since the geometry constraint (dominated by 0) is present 
at all phase angles the conclusion that high-a measure- 
ments are more effective at identifying contact binary 
systems holds. 

Figure [Hd shows how the lightcurve shape depends on 
the relative size and separation of the binary compo- 
nents. Four systems were considered - their properties 
are detailed in Table 1. Except for system #3, meant 
to represent a system formed according to the model by 
Cuk mentioned in §1, all systems are Roche binaries in 
hydrostatic equillibrium. Systems #3 and #4 produce 
shallower lightcurves due to their asymmetric mass ra- 
tio, and although the maxima and minima take up a 
large fraction of the rotational phase, the probability of 
measuring the total range of variation is not as high as 
for the more symmetric binaries #1 and #2. In conclu- 
sion, high-a observations will more easily detect contact 
binaries with similar sized components. 

Currently known binary NEOs are generally asym- 
metric with p rimary to secondary size ratios of 2 to 4 
(|Pravec et al. 2006). Furthermore, no very close pairs 
or contact binaries have been directly observed [with 
the p ossible exception of 25143 Itokawa (^Scheeres et al. 
IIOOT)]. The minimum orbital period is P = 11 hr. 
Whether these are intrinsic features or the result of ob- 
servational bias is still unclear, but the aptness of the 
method presented here to find symmetric contact bina- 
ries should help throw light on the subject. The intrinsic 
fraction of close and contact binaries is a powerful con- 
straint for models of NEO binary formation and destruc- 
tion. 

2.2. Lightcurve Phase Shift and Spin Direction 

Increasing the phase angle a produces a shift in the 
rotational phase of the lightcurves (see Fig.[2j). The shift 
happens respectively to the left or right depending on 
the object being illuminated from the left or right side, 
from the observer's standpoint. We measure a linear ro- 
tational phase shift with a slope ~ 1.4 x 10~^ per degree 
phase angle. The shift slope is similar for all four sys- 
tems in Table 1. This effect must be taken into account 
when fitting a single rotation period to data taken at dif- 
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Fig. 3. — Probability of measuring a magnitude range larger than 
Am' from random sampling each lightcurve on Fig. [2] five times. 
An observing geometry = (j) = 90° has been assumed (see text 
and Fig.[l]). 

ferent phase angles if not to be confused with evidence 
for complex rotation. 

As mentioned in §2.1, if the geometry is such that 
(j) 90° then asymmetries arise in the lightcurve. The 
direction of the asymmetry depends on the sense of rota- 
tion of the binary. Figure [U^ illustrates this effect. The 
asymmetry can thus be used to infer the sense of spin of 
the binary. 

2.3. Lightcurve Range 

Careful inspection of Fig. [2] shows that the lightcurve 
range increases with a. This is more easily seen in 
Fig. [Hi where lightcurve range is plotted versus phase 
angle (APR) for each of the four systems in Table 1. For 
comparison, a thin black solid line illustrates the same 
dependence for a triaxial Jacobi ellipsoid with axis ratios 
b/a = 0.60 and c/a = 0.43. 

We find that, with the exception of the symmetric bi- 
nary #1, the APRs of binary systems appear less reg- 
ular than that of the ellipsoid. Binary APRs show two 
roughly linear regimes, steeper for lower a and shallower 
for larger a. The slope only seems to depend on sys- 
tem type for small to intermediate phase angles: be- 
yond a ~ 50° the slopes are similar for all systems. 
Systems with smaller ranges at a = 1° have shal- 
lower ini tial APR slope s. This is similar to what was 
found by Zappala et al. ("T990) for ellipsoids. System #1, 
however, has a remarkably regular APR. The slope is 
0.004 mag/° up to a 75° and 0.001 mag/° beyond that 
value. Inciden tally, th e symmetric binary 90 Antiope 
(|Merline et all 12000: D escamps et al . 2007) has an APR 
slope 0.005 mag/ ° (Michal owski et al.ii2002t) . very close 
to that of our simulated symmetric binary. The APR is 
shape-dominated: the APR slopes do not depend signif- 
icantly on the choice of surface scattering model. 

Our results seem to indicate that less regular (more 
asymmetric) objects have steeper APRs. Our most reg- 
ular object, the symmetric binary, has the shallowest 
APR. A detailed analysis of the APR of contact binary 
systems is beyond the scope of this letter and will be 
treated elsewhere. 

3. SUMMARY 

We have presented the case for using high phase an- 
gle observations to find extreme lightcurves produced by 
contact binaries in the solar system. Sparse observations 
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Fig. 4. — a) Lightcurves of system #1 (Table 1) as a function of azimuth angle (j). Geometry is o; = 60°. b) Lightcurves of all four 
systems in Table 1 at o; = 60° and (/) = 90°. c) Asymmetric lighcurves of system #1 rotating in the prograde and retrograde directions. 
Asymmetry is due to observing geometry cf) ^ 90° and high phase angle a = 60° and shows clear dependence on the direction of rotation of 
the binary, d) Lightcurve range versus phase angle (a = 1, 5, 15, 30, 45, 60, 75, 90, 105, 120, 135, and 150°) for each system in Table 1. Thin 
solid black line corresponds to a Jacobi ellipsoid with axes ratios b/a = 0.60 and c/a = 0.43. Observing geometry is cf) = 90° (Fig. [T}. All 
four plots assume Lommel-Seeliger lunar-type surface scattering and aspect = 90°. 



of contact binaries at phase angles a > 60° are extremely 
efficient at detecting the large range of brightness vari- 
ation characteristic of these objects' lightcurves. From 
Earth, this method is most relevant to observations of 
near-Earth objects as they attain the largest phase an- 
gles. At this point, a survey of high-a measurements 
alone should probably not be used to estimate the frac- 
tion of contact binaries in the NEO population. Devia- 
tions from hydrostatic equilibrium shape and the precise 
observing geometry {0 and <p) affect the exact detection 
probability, thereby introducing uncertainty in the de- 
rived abundance. The technique is extremely efficient at 
detecting equator-on, relatively smooth contact binaries 
with symmetric components, but confirmation of the ex- 
act configuration of a system that shows large Am in just 
a few high-a measurements should always be sought us- 
ing follow-up observations. The sparse quality of planned 



all- sky surveys such as Pan- STARRS implies that a full 
shape solution can be obtained only after several years of 
data have been collected. It is therefore extremely useful 
to develop fast techniques that allow the identification 
of potentially interesting targets from a Pan-STARRS- 
type survey (or even from dedicated surveys using sub- 
Im telescopes), which can be followed up on other tele- 
scopes. 
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